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Abstract 

Nucleon matrix elements of various four-quark operators are evaluated in quenched 
lattice QCD using Wilson fermions. Some of these operators give rise to twist-four 
contributions to nucleon structure functions. Furthermore, they bear valuable infor- 
mation about the diquark structure of the nucleon. Mixing with lower-dimensional 
operators is avoided by considering appropriate representations of the flavour group. 
We find that for a certain flavour combination of baryon structure functions, twist- 
four contributions are very small. This suggests that twist- four effects for the nucleon 
might be much smaller than m'i/Q'^ . 



1 Introduction 



The knowledge of higher quark and gluon correlators in hadrons is of funda- 
mental interest in order to understand the structure of baryons and mesons on 
the basis of QCD. Matrix elements of four-quark operators contain informa- 
tion on the quark and diquark structure of the nucleon. Within the operator 
product expansion (OPE), four-quark operators give rise to higher- twist con- 
tributions (cat's-ear diagram). While this has been known for many years 
[1-6], the size of these contributions is still uncertain. Because the structure 
function ^2(0;, Q"^) of the proton is one of the best measured hadronic quanti- 
ties, the natural idea would be to extract the higher-twist contribution from 
the dependence of F2. This has, however, proven to be a difficult task (for 
a recent attempt see Ref. [7]). 
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A computation of higher-twist effects from first principles is possible by means 
of Monte Carlo simulations of lattice QCD, and first estimates using this 
method in the case of the pion became available recently [8] . In this paper we 
shall extend our previous work to the case of the nucleon. 

By means of the OPE F2{x,Q^) can be expressed through forward nucleon 
matrix elements of local operators. In the deep inelastic limit ^ oo it is 
dominated by the leading twist-two contributions. These have been the sub- 
ject of intensive studies in the past. The next-to-leading contributions have 
twist four and are suppressed by one power of More precisely, the 

OPE relates (Nachtmann) moments [9] fn dxx^^^Fofx, Q^) , which 

^ ^ L J JU ^\ ,^ J Nachtmann' 

take into account the effects of the finite proton mass nip, to the product of 
Wilson coefficients and hadronic matrix elements. Schematically one finds for 
n = 2,4,6,... 



Nachtmann 



r dxx''-^F2(x,Q^ 
Jo 

= ^C^^F2{x, Q^) [n^ + 2n + 3 + 3(n + + Amy/QY'' 



+ n{n + 2)4m;xyQ^\/{n + 2)(n + 3) 
with i = 2x/{l + ^l + 4m2a;Vg2). 

The reduced matrix elements A^*^ of twist t and spin n depend on the renor- 
malisation scale fi. The mass dimension of A^^ is t — 2. The dimensionless 
Wilson coefficients c^^ can be calculated in perturbation theory. In the fiavour- 
nonsinglet channel, the twist- two operators are two-quark operators, 

H^,D^,.-.D^'il;, (2) 
symmetrised in all indices and with trace terms subtracted. 

The four-quark operators we are interested in have twist four and higher. In 
particular, the twist-four, spin-two matrix element is given by (indices in 
{. . . } are symmetrised) 

I ^(P, S\Al^y - trace|P, S) = {PlA^^y - trace|P) = 2A^^\P^P, - trace) 
^ s 

(3) 

with the four-quark operator 

= (VSG7^75i» ii^G^.l.m (4) 
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(using the nomenclature introduced in Ref. [8]). The quark field i/j carries 
flavour, colour, and Dirac indices, the matrices t"" are the usual generators of 
colour SU(3)c, and for two flavours the flavour matrix G reads 

G = diag(e„, ea) = diag(2/3, -1/3) (5) 

in terms of the quark charges Cg. The proton states with momentum P and 
spin vector S are normalised such that 

(P, S\P', S') = {27rf2EpSiP - P')5ss' (6) 

The Wilson coefficient reads [4,5] c^^^ ^ g'^ {1 + 0{g^)). 

These expressions are to be compared with their twist-two counterparts: 

{P\0{^u} - trace|P) = 2Af\p^P, - trace) (7) 
with the operator 

= ^^G\D,ilj (8) 
and the Wilson coefficient c^2^ — 1 + 0{g^). 

The operators (4) and (8) transform identically under Lorentz transforma- 
tions, but (4) has dimension six, whereas (8) has only dimension four: four- 
quark operators will in general mix with two- quark operators of lower dimen- 
sion. This fact complicates the investigation of four-quark operators, because 
the mixing with lower- dimensional operators cannot be calculated reliably 
within perturbation theory. A nonperturbative computation in lattice QCD 
could proceed along the same lines as in the case of the twist-three matrix 
element d2 [10]. For the time being, we do not attempt such a nonperturba- 
tive calculation of the renormalisation and mixing coefficients of four-quark 
operators. Instead we restrict ourselves to cases where mixing with lower- 
dimensional operators is prohibited by flavour symmetry. 

We shall present results obtained in the quenched approximation of lattice 
QCD with Wilson fermions. A preliminary account of some of these results 
has already been given in Ref. [11]. Since the lattice formulation of gauge 
theories requires an analytic continuation from Minkowski space to Euclidean 
space, we now switch to the Euclidean formulation. In particular, all operators 
will be written down in Euclidean space, unless otherwise noted. 
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2 The general framework 



In our previous publication [8] we have studied four-quark operators in the 
pion. In this case we could avoid mixing with lower-dimensional operators by 
working with operators which carry isospin 1 = 2. Obviously, operators with 
1 = 2 vanish in the proton. Therefore wc enlarge the flavour symmetry group 
from SU(2)f to SU(3)f assuming three quarks of the same mass. Correspond- 
ingly, we then have 

G = diag(e„, ea, e,) - diag(2/3, -1/3, -1/3) (9) 

and the flavour structure of the operator in the OPE is now 

O = {tuuu + Cddd + esss){euuu + eadd + CsSs) . (10) 

While two-quark operators transform under SU(3)f according to 3®3 = 1®8, 
wc have for four-quark operators: (3(g)3) (g) (3(g)3) = 2- 1©4-8©10©T0®27. 
Four-quark operators with / = 0, 1, Is = 0, and hypercharge Y = belonging 
to the multiplets 10, 10, 27 do not mix with two-quark operators and do not 
automatically vanish in a proton expectation value. The operators belonging 
to the 27 multiplet are (giving only the flavour structure) 

Oif =^(e^ -el- 2e„e, + 2edes)[iuu){uu) - {dd){dd) 

- {us){su) - {su){us) + {ds){sd) + {sd){ds) 

— {ss){uu) — {uu){ss) + {ss){dd) + {dd){ss)] , 



^27° =F7T(e^ + + e^e^ - Se^e^ - 36^6^ + 3e^) 

[2{uu){uu) + 2{dd){dd) + {dd){uu) + {du){ud) 

+ {ud){du) + {uu){dd) — 3{us){su) — 3{su){us) (''"^•^ 

— 3{ds){sd) — 3{sd){ds) — 3{ss){uu) — 3{uu){ss) 

- 3{ss){dd) - 3{dd){ss) + 6{ss){ss)] . 

Inserting the values of the quark charges one flnds 

el -el- 2eues + 26^6^ = e^e^ - 3e„es - 36^6^ 3e^ = 1 . (13) 

As the operators belong to the same multiplet, the Wigner-Eckart theorem 
tells us that the proton matrix elements of these two operators are proportional 
to each other: 

(P|0nP) = ^(P|Ori^)- (14) 
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Fig. 1. Plateau for the bare matrix element of C^t''^ operator V^^ — trace (see 

Eq. (20)), divided by m^, with k = 0.1515. 



Furthermore, the Wigner-Eckart theorem relates proton matrix elements to 
neutron matrix elements. Thus our results can easily be rephrased in terms of 
neutron expectation values. However, unless otherwise stated, we shall only 
present the proton results. 

The operators belonging to the multiplets 10 and 10 read 

0{q^ —{dd)(uu) — {du)(ud) + {ud){du) — {uu){dd) 

+ {uu){ss) + {ds){sd) - {dd){ss) - (us)(su) (15) 
— {sd){ds) + {su){us) + {ss){dd) — {ss){uu) , 



={dd){uu) + (du){ud) - (ud){du) - {uu){dd) 

+ {uu){ss) - {ds){sd) - {dd){ss) + {us){su) (16) 
-f {sd){ds) — {su){us) + {ss){dd) — {ss){uu) . 

Being antisymmetric with respect to the interchange of the two quark-anti- 
quark pairs they do not appear in the flavour decomposition of the OPE 
operator (10). 

In the following we present Monte Carlo data from quenched simulations at 
f3 = 6.0 with Wilson fermions on a 16^ x 32 lattice. We have performed simu- 
lations at three different values of the hopping parameter k = 0.1515, 0.1530 
and 0.1550 and we have collected about 300 configurations. The statistical 
errors have been determined by the jackknife method. 
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The proton matrix elements are computed in the standard fashion from ratios 
of three-point functions, {B{t)O{T)B{0)) , to two-point functions, {B{t)B{0)), 
with the interpolating fields B and B of Ref. [12]. For <C r <C t this ratio 
should be independent of t: 

_ {B{t)0{T)Bm _ 1 

^ - {B{t)Bm - + • • • 

If we vary t keeping r fixed we should therefore find a region where R is 
constant, i.e. shows a plateau. An example of such a plateau is shown in 
Fig. 1. We have always chosen r = 5 and the spatial components of the proton 
momentum have been set to zero. The ratio R equals the matrix element in 
the lattice normalisation of states and fields. In order to obtain the fields in 
the continuum normalisation we have to multiply each quark field by ■\/2k. To 
normalise the states according to Eq. (6) we must multiply R by an additional 
factor of 2mp. We determine the matrix elements by fitting the ratio i? to a 
constant in the region 11 < i < 17. 

For a general four-quark operator the three-point function {B{t)O{T)B{0)) 
consists of three types of contributions, which can be represented pictorially 
by the following diagrams 




It is precisely through contributions of the form of the first two diagrams that 
the mixing with lower-dimensional operators occurs. Therefore these contri- 
butions cancel in the operators which we consider, and we are left with the 
contributions of the last type only. For proton matrix elements only some 
terms of the operators contribute, e.g. the {dd)(dd) terms vanish as those con- 
taining s quarks. Therefore the expectation values of the operators (11), (15), 
(16) reduce to 

{P\Oif\P) = ^{P\{uu){uu)\P)\"\ 

{P, S\Oif\P, S) = (P, S\{dd){uu) - {du){ud) + {ud){du) - {uu){dd)\P, S)\^° , 

{P,S\O^^P,S) = {P,S\{dd){uu) + {du){ud) - {ud){du) - {uu){dd)\P,S)\^'^ . 

(18) 
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1/K^ 1/K 

Fig. 2. Chiral extrapolation of the bare proton matrix element of £'27''^ for the 
operator V44 — trace (see Eq. (20)), divided by rUp. 



At each k value the matrix elements are made dimensionless by dividing 
them by the corresponding value of m^. As the bare quark mass is given 
by aniq = 1/ {2k,) — 1/ (2kc) we perform the chiral limit by linear extrapola- 
tion in l/fi;to 1/k, = l/nc = 6.3642. As an example of our results we show in 
Fig. 2 the chiral extrapolation of the bare proton matrix element of 1^44 — trace 
(7 = 1 component in the 27 representation of SU(3)f) divided by m^. (For the 
definition of V^fj, see Eq. (20).) 



3 Operators from the 27 multiplet 

The twist-four contribution in the F2 structure function comes from the four- 
quark operator A^^, see Eq. (4). In order to access the flavour-27 component 
experimentally one has to combine the structure functions of several baryons 
(p, n, A, E, S) in such a way as to project out the desired flavour combination, 
e.g. 

{j)\Oif\p) = (S+|C|S+) - 2(E°|0|S°) + 
= -(S+|C>|E+) - 
-6(A|C|A) + 2(S°|C|S°) + 2{E-\0\E-) 
+2{p\0\p) + 2{n\0\n) . (19) 
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Table 1 

Matrix elements of the spin- two operators from the 27 multiplet, divided by 
and extrapolated to the chiral limit. 



operator 





— trace 


(-0.9 ± 0.8) • 10-4 




— trace 


(-4.5 ±1.0) • 10-4 




— trace 


(4.9 ±0.8) • 10-4 




— trace 


(-0.0 ± 1.3) • 10-4 




— trace 


(8.4 ± 1.4) • 10-4 




— trace 


(-4.6 ±1.9) • 10-4 



Unfortunately most of these terms will not be measured in the foreseeable 
future. So a direct comparison with data is out of question. On the other 
hand, they can be used as a testing ground for models of hadrons, taking the 
role of experimental data. Note that the 27 contribution can also be isolated 
by studying combinations of electromagnetic and weak structure functions 
[13]. 

Of course, we need to know the renormalised operators. Although due to our 
choice of the flavour-27 component mixing with two-quark operators is absent, 
different four-quark operators may still mix under renormalisation. Therefore 
we have computed the matrix elements of the following operators (using the 
nomenclature introduced in Ref. [8]): 



' fJ.U 








rpc 










= '0<^7/.75^ V'G'7i.75^, 


T 


= ijjGanp'i/j'ijjGai^pijj. 



The bare expectation values divided by rrip and extrapolated to the chiral 
limit are given in Table 1 for the the spin-two components, while the traces 
are given in Table 2. E.g. the number shown for the operator A^^^ — trace in 
Table 1 is what we obtain for (P|^('U7475t"M)('U7475t"M) — trace|P)/m4 in the 
chiral hmit. We have checked that these operators fulfil their Fierz identities. 
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Table 2 

Matrix elements of the spin-zero operators from the 27 multiplet, divided by 
and extrapolated to the chiral limit. 



Dirac structure 




ta^ta 


1(8)1 


75 ® 75 


(2 


8 ± 0.7) • lO-'' 


(-13.4 ±2.6) • IQ-^ 


10 1 


(- 


7.1 ± 1.2) • IQ-^ 


(17.8 ± 2.9) • IQ-'' 


7/^75 <^ 7m75 


(17.2 ±2.0) • 10-^ 


(-12.1 ±3.9) • 10"^ 


7/*®7m 


(- 


18.2 ±2.2) • 10""^ 


(7.4 ±4.0) • 10-"^ 




(- 


4 ± 7) • 10"^ 


(17.6 ± 9.6) • 10"^ 



The renormalisation constants have been calculated in one-loop perturbation 
theory [8]. The renormalised spin- two piece of the operator reads 



9i 



KM =Al.-T^ (3 ln(a/x) +46.072285) A 



167r2 



/ 8. , . 

+ 



- ln(a//) + 0.083982 j V^^, ^^^^ 

+ (-^ln(a/x) + 0.157467) 
- 1.071448 T^^ - 2.008965 T^^ 

where go is the bare coupling constant (/9 = 6/g'Q).The renormalisation scale 
H will be identified with the inverse lattice spacing 1/a. In our simulations 
this has a value of 1/a ~ 2.12 GcV (using ro to set the scale). In terms of the 
renormalised operator the reduced matrix element 742^'* is then given by 



1 (4) 27,1=1 _ 2 (-P|^( ^7475^°^)(^7475^°^) -trace|P) 

9 2 

mi 



3 m4 



(22) 



and we obtain for the lowest moment of F2 in our special flavour channel 

/' dx F^ix, Q') = -0.0005(5) ^^^4^ + O(a^) . (23) 

Jo ^ Nachtmann ^ ' Q2 V V / 

The analogous result for the neutron differs from the above only by the sign. 

In the proton the corresponding twist- two contribution is about 0.14 at — 
5 GeV^. As in the pion, the twist-four correction is tiny. Our result may be 
compared with bag model estimates. In this model the scale for the prefactor in 
Eq. (23) is set by B/m^ ^ 0.0006, where B ^ (MSMeV)"^ is the bag constant. 
The factor B/rrip is however multiplied by a relatively large (and negative) 
number [3]. 
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It is rather difficult to determine the first moment of the higher-twist contri- 
bution to ^2(2;) experimentally. Phenomenological fits to the available data 
give a positive value of about 0.005(4) GeV^/Q^ [^,14]. Our matrix element, 
which is due to its flavour structure only one contribution to the full moment, 
is considerably smaller than this phenomenological number. 



4 Operators from the 10 and 10 multiplets 

Having found rather small matrix elements for our four-quark operators from 
the 27 one may ask if operators from the 10 or 10 of SU(3)f (although not 
contributing to F2 in the OPE) would have larger matrix elements. With the 
two possible colour structures that can form colour singlet operators, these 
operators are linear combinations of terms of the form {ipG rt''ilj){ijjG' T't'^'tp) 
and {'4^GT'4>){'4>G'V''4>), respectively, where F and V are Dirac matrices. We 
have chosen the flavour matrices G, G' such that we get the following flavour 
structures: 

(dd)(uu) - (uu)(dd) (24) 

and 

{du){ud) - {ud){du) . (25) 

These can be combined to yield the 10 and 10 structures in Eq. (18). 

Discrete symmetries impose restrictions on the matrix elements of these op- 
erators. We have 

(P, S\0\P, Sy = (P, -S\TVOV-^T-^\P, -S) = (P, S\O^P, S) (26) 

where V is the parity and T is the time inversion operator. For the Dirac 
matrices used in our computations we define sign factors si, s'l, S2 and s'2 by 

74Ft74 = siF , 74F't74 = s;f' , ^^^^ 

7475CF*C-i7574 = S2T , 7475CF'*C-i7574 = s'^V . 

Here G is the charge conjugation matrix with G^'^G~^ — —^^i- One more sign 
eo is determined by 

(P, -S\0\P, -S) = eo(P, S\0\P, S) . (28) 
Prom Eq. (26) we now get for the fiavour structure (24) 

(P, S\0\P, sy = eoS2s'2{P, S\0\P, S) = sis[{P, S\0\P, S) (29) 
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Table 3 

Expectation values of operators with the flavour structures (24) and (25) in an 
unpolarised proton, divided by nip and extrapolated to the chiral limit. 



Dirac structure 


flavour 


taSta 


1(8)1 


1 (g) 74 


(24) 


(0.4 ± 0.6) • 10-3 


(3.0 ± 1.6) • 10-3 


e4o/357a75 'Si a ps 


(24) 


(-5.6 ±2.8) • IQ-^i 


(5.9 ±3.6) • lO-^i 


75 ® 7475 


(25) 


(0.6 ±0.5) • W-^ 


(-3.3 ± 1.8) • 10-3 


7a <^ 0-4a 


(25) 


(-3.1 ± 1.2) • IQ-^i 


(0.1 ± 1.9) • 10-3i 



and for the flavour structure (25) 

{P,S\0\P,S)* = eoS2s'2{P,S\0\P,S)^-sis[{P,S\0\P,S) . (30) 

Thus the matrix element {P, S\0\P, S) is real if €08282 — 1 and purely imag- 
inary if €08282 = —1; the matrix element vanishes if €08282 — — Sis^ for the 
flavour structure (24) or €08282 = sis^ for the flavour structure (25). We have 
checked that these restrictions are satisfied by our results within statistical 
errors. We restrict ourselves in the following to the matrix elements which are 
not forced to be zero by the above relations. Note that for a given Dirac struc- 
ture at most one of the flavour structures (24) and (25) yields a non-vanishing 
result. 

The definite Lorentz transformation properties of our operators could be used 
to define reduced matrix elements, e.g. in Minkowski space one gets 

(P, S\{d^^^,d){u^,j,u) - iu^^j,u)idj,^,d)\P, S) = Ae^,a^(P">S^ - S'^P") . 

(31) 

Thus in this case the matrix element with /j, = 1, 1/ — 2 and 5"^ = Sa3 is equal 
to the one with /jl — 2, u — 3 and S"' — Sai- This holds only on average, so 
in order to increase the statistics we averaged over these matrix elements to 
reduce the statistical error. The bare expectation values divided by are 
given together with their statistical errors in Tables 3 and 4. 

The order of magnitude of the results does not differ greatly from those found 
for the operators in the 27. The renormalisation constants for the 10 and 10 
operators are not known, but we do not expect that the rcnormalised operators 
have much larger matrix elements than the bare ones. 
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Table 4 

Expectation values of operators with the flavour structures (24) and (25) in a po- 
larised proton (S = 63), divided by and extrapolated to the chiral limit. 



Dirac structure 


flavour 


ta^ta 


101 


1 <8) 7375 


(24) 


(-3.3 ±0.6) • 10-^i 


(15.9 ± 1.6) • 10-3i 


1 (721 


(24) 


(-4.1 ±0.6) • 10-3 


(10.7 ± 1.8) • 10-3 


75 ® 0"43 


(24) 


(2.2 ±0.5) • 10-3 


(-4.4 ± 1.3) • 10-3 


74 ® 7375 - 73 <^ 7475 


(24) 


(-9.1 ±0.9) • lO-^i 


(13.0 ±2.1) • 10-3i 




(24) 


(42 ± 4) • 10-3 


(-51 ± 5) • 10-3 


75 <^73 


(25) 


(4.9 ± 0.5) • 10-3i 


(-15.6 ± 1.4) • 10-3i 


72 ® 71 


(25) 


(2.4 ± 0.5) • 10-3i 


(-5.6 ± 0.9) ■ 10-3i 


7275 ® 7175 


(25) 


(6.1 ± 0.8) • 10-3i 


(-7.6 ± 1.1) • 10-3i 


7q75 <^ crsa 


(25) 


(21.5 ± 1.8) • 10-3 


(-19.2 ±2.6) • 10-3 


(72a ® O'la 


(25) 


(-8.9 ± 1.5) • 10-3i 


(8.1 ± 1.8) • 10-3i 



5 Diquarks 

The four-quark operators can be rewritten to look like a diquark density. We 
have computed matrix elements of operators of the following form: 

^{UaI'l5Cli^){u^,C~^-f5^'ub'){5ab'Sba' " Saa'Sbb') , (32) 
Y^{Ua^l5CuJ){u^,C~^-f5r'ub'){Sab'Sba' + ^aa'hv) , (33) 

where a, b, a' and b' are the colour indices. These are the two possibilities 
to form a colour singlet. In Eq. (32) the diquark is in a 3 of colour and thus 
anti-symmetric in colour. In Eq. (33) it is in a 6 and symmetric in colour. 
Because of the Pauli principle it has to be symmetric (anti-symmetric) in the 
other indices. The flavour structure being symmetric, the Dirac structure has 
therefore to be symmetric (anti-symmetric). Thus a given Dirac structure will 
appear only for one of the two possible colour structures. The expectation 
values of operators of the form (32) and (33) can be computed from those of 
the four-quark operators studied in Section 3. But in order to get the correct 
errors we have redone the analysis. The results (again divided by m^) are 
presented in Tables 5 and 6 for the spin-zero and the spin-two contributions, 
respectively. 



12 



Table 5 

Matrix elements of the spin-zero operators from the 27 multiplet in the diquark 
picture, divided by nip and extrapolated to the chiral limit. 



operator 


colour 




Jo {u'yf,'y5Cu^){u^C''^-f5-f^u) 


3 
3 
6 
6 
6 


(-83 ± 8) • lO-'^ 
(42 ± 22) • 10"^ 
(45 ± 9) • 10-'^ 
(-0.8 ± 2.6) • lO"-^ 
(-4.9 ±6.7) • 10-"^ 



Table 6 

Matrix elements of the spin-two operators from the 27 multiplet in the diquark 
picture, divided by mf, and extrapolated to the chiral limit. 



operator 


colour 




3^ (12747575 Cii'^)(tt'^C"^757475u) - trace 
jo{ua4a75Cu^){v7- C^^^^a^au) — trace 
(117475 C'u'^)('u'^C"^7574u) - trace 


3 
3 
6 


(12.9 ±2.2) • 10-^ 
(-16.8 ±4.3) • 10-"^ 
(3.9 ±2.9) • 10"'' 



Strictly speaking, we are again studying operators from the 27 representation 
of SU(3)f, whose uuuu component is given by Eqs. (32) and (33), respectively. 
At least within the quenched approximation it seems however reasonable to 
consider (32) and (33) as representing valence diquark densities. In the same 
spirit, one could also investigate ud diquarks. But due to flavour symmetry 
(see Eq. (14)) the corresponding matrix elements are proportional to those of 
the uu diquarks (32) and (33). Writing down only the flavour structure one 
flnds for the expectation values in the proton 

{P\{ud'^){d^u)\P) = \{P\{uu^){u''u)\P). (34) 

In order to interpret our results we have combined operators from Tables 5 and 
6 such that they correspond to diquarks of spin zero and spin one. Speciflcally, 

for an operator O^j, with two space-time indices wc take the expectation value 
of to represent a spin-zero diquark and the expectation value of J2i=i 
to correspond to a spin-one diquark. The results (once again completely re- 
analysed) are given in Table 7. 
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Table 7 

Expectation values of operators from the 27 multiplet which correspond to diquarks 
of spin zero and spin one, divided by and extrapolated to the chiral limit. For 
the operator O^^, the spin-zero contribution is O44, the spin-one contribution is 



operator 


colour 


spin 




spin 1 




3 
3 
6 


(-6.8 ± 2.4) 
(-5.7 ±6.1) 
(15.0 ±3.6) • 


• 10-^ 

• lo-'' 

IQ-^ 


(-76 ± 7) • IQ-^ 
(48 ± 18) • lO-'' 
(30 ± 7) • 10-^ 



For the operators in the 3 of colour the absolute values for the spin-one di- 
quarks are considerably larger than for the spin-zero diquarks. For the single 
operator in the 6 of colour the difference is less pronounced. This pattern 
can tentatively be understood in a non-relativistic quark picture. When the 
diquark is in the 3 of SU(3)c it is anti-symmetric in the colour indices, and 
therefore the symmetric (in the spin indices) spin-one state is favoured over 
the anti-symmetric spin-one state. On the other hand, when the diquark is in 
the (symmetric) 6 of colour one might at first sight expect the anti-symmetric 
spin-zero state to dominate over the symmetric spin-one state. Although the 
spin-zero contribution is indeed less suppressed than in the 3 case it is not 
really dominating. This is probably related to the fact that a diquark in the 6 
of SU(3)c must be accompanied by (at least) one gluon if it is to form a colour 
singlet with the remaining quark. The coupling to the gluon, mixing "large" 
and "small" components of the quark spinors, would invalidate the above ar- 
guments which worked reasonably well for diquarks in the 3 of colour. 

Of course, the operators from the 10 and 10 multiplets can also rewritten in 
diquark form. They then appear as linear combinations of operators in which 
the diquark is either in a 3 of SU(3)c or in a 6. The fact that the matrix 
elements in Table 4 for the colour structures ta ta and 1 (g) 1 have opposite 
signs translates into a suppression of the 6 diquarks relative to the 3 diquarks. 
This is in accord with our observations made above in the case of the operators 
from the 27 multiplet. 



6 Summciry 

In this paper we have computed the expectation values of a variety of four- 
quark operators in the proton by means of quenched Monte Carlo simulations. 
Since it is rather difficult to treat the mixing with lower-dimensional operators 
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correctly, wc have restricted ourselves to operators whose flavour structure 
forbids this kind of mixing. The additional requirement that the operators 
should not automatically vanish in the nucleon led us to consider the flavour 
group SU(3)f and to choose operators from the 27, 10, and 10 representations. 
One of the operators from the 27 is responsible for the twist-four contribution 
to the lowest moment of F2 in a somewhat exotic flavour channel. We find 
a rather small value for this contribution. This observation disfavours large 
higher-twist effects in general, although we cannot exclude, of course, that our 
result is due to strong cancellations between different flavour contributions. 

Still, and with all due caution, our flndings flt into a general trend emerging 
from various pieces of information on higher-twist contributions to moments 
of nucleon structure functions: The natural energy scale for the correspond- 
ing correlators lies well below the nucleon mass leading to small numerical 
coefficients when the matrix elements are expressed as multiples of (mp)*~^. 

Thus we arrive for the nucleon at a conclusion which is similar to what we 
observed in the pion [8] . For a more detailed comparison we plot in Fig. 3 the 
renormalised pion matrix elements [8] with the fiavour structure 

{t:^\{uu){uu) + {dd){dd) — {uu){dd) — {M){uu) 

- {ud) {du) - {du) {ud) |7r+) /ml ^ ' 

together with the corresponding renormalised matrix elements for the proton 
(p|10 ■ 0^^^\p)/'m'^ (in lattice units). We display the results for the spin-two 
components setting = i/ = 4 (with the trace term subtracted). The nor- 
malisation of the operators is chosen such that the fiavour structure {uu){uu) 
appears with the factor 1 in both cases. (Alternatively, it may be remarked 
that SU(3)f makes the above pion matrix element equal to the expectation 
value of 10 ■ in the meson-octet analogue of the proton, the .) It is 

no great surprise that the numbers do not show many similarities - after all, 
the pion and the proton are very different particles. 

Four-quark operators in the 10 and 10 representations do not lead to much 
larger matrix elements than the 27 operators, although quite a few of those 
which we studied give clean signals. In the 27 sector, a rewriting of our oper- 
ators in terms of diquarks reveals a structure which lends itself to an interpre- 
tation with the help of quark model ideas: Diquarks in the 3 of SU(3)c have 
preferably spin one. Diquarks in the 6 of SU(3)c, on the other hand, do not 
fit so well into a nonrelativistic picture. 

Higher-twist effects will challenge lattice QCD for a few more years. Our in- 
vestigations show that four-quark operators can give reasonable signals in 
present quenched Monte Carlo simulations. But the study of physically more 
interesting flavour channels and further twist-four operators remains an open 
problem whose solution requires progress in nonperturbative renormalisation. 
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Fig. 3. Renormalised four-quark matrix elements in the pion and in the proton (in 
lattice units). 



especially in the treatment of mixing with operators of lower dimension. 
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